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ABSTRACT 


The effects of ionizing radiation on selected liquid propellants have 
been investigated and preliminary development of a comprehensive computer pro- 
gram for shield optimization has been completed. Gamma -ray attenuation param- 
eters for six candidate propellant materials have been calculated. These data 
are formulated into linear attenuation coefficients and buildup factors which 
are tabulated. A logic flow diagram for generalized weight optimization to 
select shield configurations for various propellant materials was generated 
and a simplified computer program which minimizes the weight of a primary 
direct laminated shield was developed. Radiation effects on three fuels and 
three oxidizers were investigated in several container materials and over a 
range of dose rates. Propellant materials investigated include; FLOX, oxygen 
difluoride, nitrogen tetroxide (inhibited and uninhibited), diborane, hydra- 
zine, and the LPG fuel, propane. In general, it was found that oxidizers are 
much less susceptible to radiation induced decomposition than fuels and 
differences due to dose rate are generally not large. 


RADIATION EFFECTS ON LIQUID PROPELLANTS 


INTRODUCTION 


It is advantageous to use high specific impulse liquid propellants 
on long-term space probes using isotope or nuclear power sources for auxiliary 
power because of the rigid weight and volume requirements on such missions. 
However, more definitive knowledge of environmental effects on the long-term 
storability of these propellants is desirable so that design engineers may 
have reasonable assurance that the stored propellant system chosen for a 
particular mission will fulfill its function successfully. A preliminary 
study under Contract NAS7-577 was carried out to ascertain the state of the 
art regarding radiation effects data on liquid propellants and to analyze the 
suitability of using various propellant systems on space probes with different 
auxiliary power sources. 

The present program was established with the objective of developing 
more detailed Information on radiation effects and shielding requirements for 
a few selected liquid propellants. Three tasks were established in this pro- 
gram. Task 1 was to obtain radiation effects data on three fuels and three 
oxidizers irradiated as liquids. Task 2 was to begin development of a com- 
puter program to ascertain what special design criteria or shielding tech- 
niques must be employed to protect components of the propulsion system from 
radiation, if necessary. Task 3 was a Phase 1 storage test in a low dose 
rate radiation field extending over a period of months. Results on all three 
tasks are detailed in the following sections. Task 2, Radiation Shielding 
Optimization, is discussed in the first section of the report. Tasks 1 and 3 
are then discussed in the following sections. Results from both Tasks 1 and 3 
are reported together for a given material so that results may be more easily 
compared and evaluated. Wherever possible, data have been presented in 
engineering units to facilitate their use among the widest possible audience. 
Where units specific to radiation effects studies, such as G-values, have been 
used, an explanation of the term and its use has been given. 


RADIATION SHIELDING OPTIMIZATION 


Task 2 

In a study previously conducted by Battelle's Columbus Laboratories 
for the National Aeronautics and Space Administration (NAS7-577), the probable 
types, spectra, and intensities of radiation from nuclear power sources which 
may possibly be used in deep space missions were determined. The approach 
followed in the study was a conservative one, i.e., that of determining radia- 
tion intensities on the basis of an extrapolated point source without consid- 
eration of the self-shielding inherent in a nuclear power source of finite 
geometry. The results of the study indicated that radiation damage to liquid 
propellants from certain types of isotopic heat sources, viz., ^^°Pu 02 and 
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^'^Cm 203 , which are primarily alpha emitters, would most likely* be insignifi- 
cant even in unshielded cases for source power levels up to 200 kw (thermal) 
and corresponding mission lifetimes extending to 5 years. Conversely, the same 
study revealed that in the case of other power sources investigated, e.g., 
90srTiO3, ^*^€0 (metal) and low-powered nuclear reactor sources such as the 
SNAP-lOA reactor, radiation damage to liquid propellants could be a signifi- 
cant problem and, therefore, should be further investigated from the standpoint 
of reducing the radiation environment. It must be remembered, however, that 
the significance of a given radiation effect is largely determined by the pro- 
posed application of the material under investigation. Thus, as described in 
later sections of this report, the changes in chemical composition resulting 
from irradiation of the propellant materials investigated were in general very 
small and probably not significant from the standpoint of their chemical prop- 
erties. However, in most cases such small chemical changes were accompanied 
by the accumulation of noncondensible gaseous products of the liquid propellant 
samples. The continual buildup of pressure over the liquid propellant could 
present severe engineering problems, particularly on long missions, even 
though the total liquid decomposed or chemically altered was small. 

Attenuation of the radiation intensity in the vicinity of the pro- 
pellant can be accomplished by (1) the use of a radiation shield and/or 
(2) an increased separation distance between the nuclear power source and the 
liquid propellant. However, overall system requirements (i.e., weight and 
configuration restrictions, etc.) may limit the extent to which either of the 
foregoing techniques may be utilized; therefore, a method of systems analysis 
is required which will consider the nuclear power source, shield, and pro- 
pellant subsystems within the system restraints while insuring that radiation 
damage to the propellant does not reach some critical level. 

The ultimate selection of optimum liquid propellants for deep space 
missions may well be predicted on the shielding properties of these materials. 
Consider, for example, the case where a particular propellant is initially 
chosen on the basis that it requires only a minimum shield in the close prox- 
imity of the nuclear power source while receiving an integrated dose below its 
radiation damage threshold. Further analysis, however, of the overall system 
restraints may well indicate the necessity to increase the initially determined 
shield thickness and/or the source-to-propellant separation distance. A re- 
straint of the latter type would be the allowable integrated dose for the 
transistors in the instrument package (subsystem) which might possibly be 
exceeded for the shield and geometric configuration determined above. Therefore, 
an increase in either of the two parameters, viz., shield thickness and source- 
to-propellant separation distance, in an effort to preclude the possibility of 
exceeding the dose limitation for transistors will correspondingly reduce the 
integrated dose encountered by the propellant, and may possibly permit the 
selection of a liquid propellant which had previously been eliminated on the 
basis of radiation damage. 


* Due to a lack of experimental data for the neutron irradiation of liquid pro- 
pellants, there is some uncertainty concerning the radiation effects resulting 
from (n, a) reactions . 
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On the basis of the foregoing considerations, an optimum computer 
program should contain a library* ** of neutron and gamma radiation attenuation 
data whereby it would account for the effects of auxiliary (propellant) 
shielding in the overall system optimization from the standpoint of shield 
weight . 

The results of the initial study to provide the basis for a compre- 
hensive shielding optimization computer program are detailed in the following 
sections . 


Subtask A. Compilation and/or Calculation 
of Propellant Attenuation Data 


Compilation of Photon and Neutron Cross-Section Data 

Ten energy group photon cross-section data were obtained for pro- 
pellant constituents, viz., H, N, 0, F, and C, from published sources and pre- 
pared in a format compatible with an in-house transport code, ANISN(^)*'', for 
subsequent use in propellant attenuation calculations. These data covered a 
photon energy range of 0.01 to 4 Mev and consisted of photoelectric, pair 
production, and Compton scattering cross sections; the latter were represented 
by the zeroth through third moments of the total transfer matrix for use in 
anisotropic scattering calculations. 


Calculation of Propellant Gamma Attenuation Data 

A check-case calculation was made with the ANISN code to simulate a 
benchmark experiment for gamma-ray transport through water. Water was selec- 
ted as a shield material for the first calculation because ( 1 ) experimental 
data are readily available for this material and ( 2 ) its elemental composition 
is similar to that of the liquid propellants under consideration. The input 
data utilized for the calculation included an Sg angular quadrature, Pg order 
of scattering, and spatial mesh size of 2 cm; output data were found to be 
within 10 percent of the experimental values reported for linear attenuation 
coefficient (|j,) and buildup factors. 

On the basis of the good correlation obtained between calculated and 
experimental results for water, additional ANISN calculations were made to 
determine the attenuation of gamma rays from a point source in hydrazine (N 2 H 4 ) , 
FLOX, OF 2 , LPG (CH 4 ), N 2 O 4 , and diborane (B 2 H 5 ) . Linear attenuation coeffi- 
cients and dose buildup factors for gamma-ray penetrations up to 160 cm in these 
fuels were compiled for up to seven photon energy groups (Table 1) . 


* Or else a subroutine whereby such data could be calculated. 

** References appear on page 21. 
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table 1. PHOTON ENERGY GROUP STRUCTURE 

Group No. 

Energy Range , Mev 

1 

4 3 

2 

3 - 2.5 

3 

2.5 -* 2.0 

4 

2.0 1.5 

. 5 

1.5 -* 1.0 

6 

1.0 - 0.75 

7 

0.75 - 0.5 

table 2. N 2 O 4 BUILDUP 
ATTENUATION 

FACTOR AND LINEAR 
COEFFICIENT DATA 




Penetration 

Depth 

(cm) /Buildup Factor 

- 

Group No. 


40 

80 

160 

0.054 

1 

1.2 

1.7 

2.0 

2.5 

0.061 

2 

1.3 

1.85 

2.25 

2.75 


0.069 

3 

1.4 

2.0 

2.5 

3.0 

0.077 

4 

1.45 

2.0 

2.5 

2.7 

0.085 

5 

1.5 

1.9 

2.3 

2.4 

0.105 

6 

2.0 

2.7 

3.6 

4.2 

0.114 

7 

1.9 

2.4 

2.9 

2.9 
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Sub task B. Formulation of Propellant Attenuation Data 

The dose buildup factors computed for four penetration depths (20, 
40, 80, and 160 cm) and linear attenuation coefficients computed by ANISN 
are listed by group in Tables 2 through 7. 


Subtask C. Construction of Computer 
Program Logic Diagram 


Basic Logic Diagram 

Under this task, the program logic for the radiation shielding 
optimization code was constructed. The sequential flow of the basic program 
is shown in Figure 1. Segments of the program include (1) input of problem 
data, (2) preliminary data editing and manipulation, (3) selection of cal- 
culational sequences for weight optimization (i.e., monitor control), 

(4) coarse (Level 1) optimization of shield, (5) fine (Level 2) optimization 
of shield, and (6) an output edit. Each of these program segments is described 
below. 


Data Input . Data planned for input to the computer program include 
source type and configuration, source-to-propellant-to-instrument package 
geometry, shield type (split, laminated, and scatter) and configuration, 
materials together with initial locations, propellant and instrument package 
radiation dose limitations, and program options for the calculational methods, 
shield positions, and level of optimization data edit. 


Data Edit . The data edit segment of the code will set up and store, 
in the proper format geometry, data for the various dose calculation routines, 
determine the radiation source strength and spectrum as well as arrange other 
pertinent data in a format appropriate for outlines used in the shield 
optimization. 


Control Monitor . Under control monitorship, iterations and flow con- 
trols will be set up for specific vehicle configurations. This segment of the 
program includes optimum sequencing of the various routines used in the selec- 
tion of a minimum-weight shield. 


Level 1 Optimization . The Level 1 or coarse optimization routine 
selects the weight-optimized shield configuration and position (s) through the 
use of transmission matrix tables for shield calculations. This method re- 
quires a minimum of machine (computer) time per calculation; however, the 
accuracy of the results is limited. Primary use of the Level 1 weight optimi- 
zation is to provide a starting point for the Level 2 optimization. In 
addition, the Level 1 optimization calculations can be used for parametric 
studies of space systems where a high level of accuracy is not required. 
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TABLE 3. FLOX BUILDUP FACTOR AND LINEAR 
ATTENUATION COEFFICIENT DATA 


IX 

Group No. 

Penetration Depth 

(cm)/Buildup Factor 

20 

40 

80 

160 

0.045 

1 

1.2 

1.5 

2.0 

2.5 

0.051 

2 

1.3 

1.7 

2.2 

2.7 

0.057 

3 

1.4 

1.9 

2.5 

3.0 

0.064 

4 

1.5 

1.9 

2.5 

2.7 

0.070 

5 

1.5 

1.9 

2.5 

2.5 

0.086 

6 

1.8 

2.3 

2.5 

3.0 

0.092 

7 

1.7 

2.2 

2.5 

2.9 



TABLE 4. 

0F2 BUILDUP 

attenuation 

FACTOR AND 
COEFFICIENT 

LINEAR 

DATA 




Penetration Depth (cm)/Buildu 

p Factor 

U 

Group No. 

20 

40 

80 

160 

0.048 

1 

1.2 

1.6 

2.1 

2.6 

0.054 

2 

1.3 

1.8 

2.3 

2.8 







0.061 

3 

1.4 

1.9 

2.5 

3.0 

0.068 

4 

1.5 

2.0 

2.6 

2.9 

0.075 

5 

1.5 

2.0 

2.4 

2.5 

0.093 

6 

1.9 

2.4 

2.9 

3.0 

0.10 

7 

1.8 

2.3 

2.7 

2.9 
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TABLE 5. N 2 H 4 BUILDUP FACTOR AND LINEAR 
ATTENUATION COEFFICIENT DATA 


u 

Group No. 

Penetration 

Depth 

(cm) /Bui 

Idup Factor 

20 

40 

80 

160 

0.034 

1 

1.2 

1.6 

2.2 

2.6 

0.038 

2 

1.2 

1.6 

2.2 

2.6 

0.043 

3 

1.3 

1.7 

2.3 

2.6 

0.049 

4 

1.3 

1.8 

2.4 

2.7 

0.055 

5 

1.4 

1.7 

2.0 

2.3 

0.068 

6 

1.6 

2.2 

3.0 

3.9 

0.072 

7 

1.6 

2.0 

2.5 

2.9 



TABLE 6 . 

CH 4 BUILDUP 
ATTENUATION 

FACTOR AND LINEAR 
COEFFICIENT DATA 


- - ■ ■ 


Penetration Depth (i 

cm) /Bui 

Idup Factor 

M- 

Group No. 

20 

40 

80 

160 


0,019 

1 

1.1 

1.2 

1.9 

2.5 

0.021 

2 

1.1 

1.3 

1.9 

2.6 

0.024 

3 

1.2 

1.3 

2.0 

2.6 

0.027 

4 

1.2 

1.4 

2.1 

2.7 

0.031 

5 

1.2 

1.5 

2.1 

2.8 

0.038 

6 

1.3 

1.5 

2.2 

3.0 

0.041 

7 

1.4 

1.7 

2.2 

2.7 


4c 


TABLE 7. B2H6 BUILDUP FACTOR AND LINEAR 
ATTENUATION COEFFICIENT DATA 






FIGURE 1. BASIC PROGRAM LOGIC DIAGRAM 








5 


Level 2 Optimization . The Level 2 optimization routine starts with 
the Level 1 optimization results and performs a highly accurate weight optimi- 
zation based on shielding calculations incorporating transport theory. These 
calculations will evaluate in detail the effects of flux perturbations at 
shield interfaces, the secondary radiation produced in the shield, propellant 
and structural members, and the anistropic flux distributions in these mem- 
bers. Both one- and two-dimensional transport -theory calculations will be 
required to obtain the highly accurate weight-optimized systems. 


Output Edit . In the output edit segment of the code, the calcula- 
tional data generated during the optimization are edited. All data required 
for the weight optimization are printed. 


Semidetailed Logic Diagram 

The semidetailed logic for the Level 1 and Level 2 optimization seg- 
ments of the basic computer program (see previous section) is illustrated in 
Figures 2, 3, and 4.* The diagram in Figure 2 includes a portion of the logic 
centered with shield structure (laminations) and position as well as weight 
optimization of the shield. Figure 3 illustrates the optimization logic for 
split-shield systems. Figure 4 illustrates the optimization logic for systems 
recurring scatter-shield analyses. Details of the optimization routine have 
been omitted in order to illustrate the basic logic of these segments of the 
program. 


Analytical Model on Which Program Logic is Based 

The analytical model illustrated in Figure 5 with spherical oxidizer 
and fuel tanks in tandem has served as the generalized framework for the com- 
puter program logic under construction. It is assumed that the component 
(oxidizer or fuel) having the higher resistance to radiation damage will be 
located between the remaining component and the radiation source. 


Subtask D. Preliminary Programming of Selected 
Operations of Computer Model 

The optimization of the "primary direct laminated shield" was chosen 
as the selected operation of the computer model to be programmed in this con- 
tract performance period (see enclosed portion of logic diagram in Figure 6) . 
The program is written for use in both the Level 1 and Level 2 optimization 
segments of the generalized computer program (see Figure 1) which is to be 
developed in the future. The difference between Level 1 and Level 2 uses of 
this routine will be in the transmission matrix data utilized. In the Level 1 
optimization segment, the routine will use generalized data in the form of 
point value transmission matrices to determine the transmitted particle or 


* Definitions of commonly used terms in these figures are given in Appendix A. 
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photon flux for each material traversed by the source radiation whereas, in the 
Level 2 optimization segment, the routine will rely on transmission matrices 
computed by another subroutine incorporated in the generalized optimization. 

The routine as written will determine the optimum shield configura- 
tion for the system shown in Figure 7 which represents a portion of the 
analytical model shown in Figure 5. 


Computer Program for the Optimization 
of "Primary Direct Laminated Shield" 

A listing of the computer program which utilizes the transmission 
matrix technique is presented in Appendix C of this report. A sample problem 
input and output for this program is also listed. 


SUMMARY OF TASK 2 


Gamma-ray attenuation parameters for six candidate propellant con- 
stituents have been calculated. These data were then formulated into linear 
attenuation coefficients and buildup factors which have been presented in 
tabular form. A logic flow diagram for a generalized weight optimization com- 
puter program to select shielding configurations for propellant materials was 
generated. A simplified computer program which minimizes the weight of a 
primary direct laminated shield was developed. The program uses a transmis- 
sion matrix solution technique, which reduces computation time, to calculate 
the radiation fluxes through the shield. Future work in this technical area 
should include (1) Initiation of the programming of the generalized computer 
program to optimize the weight of space-shield propellant systems and (2) form- 
ulation of transmission matrices for use with the optimization technique. 


RADIATION EFFECTS 


Tasks 1 and 3 


Experiments to obtain radiation effects data on Tasks 1 and 3 are 
described in this section of the report. Task 1 is radiation effects on 
several selected fuels and oxidizers as liquids while Task 3 is a Phase 1 
storage test in a low dose rate radiation field. The liquid fuels examined 
in both tasks are hydrazine, diborane, and an LPG fuel, propane, and the 
liquid oxidizers are FLOX, oxygen difluoride, and nitrogen tetroxide, both 
inhibited and uninhibited forms. The results obtained from the two tasks 
are reported and discussed together in each fuel material category so that 
all information on a given material is together in a single section of the 
report. 
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Throughout this section of the report results are generally given in 
terms of G-values for gas formation resulting from irradiation. Radiation 
doses are given in rads. G-values are defined as the number of molecules of 
product formed per 100 ev of energy absorbed by the material in question. The 
unit of energy absorption, the rad, is defined as 100 ergs absorbed per gram 
of sample. The data are thus presented in a very general way so that the re- 
sults can be easily applied to estimating radiation effects in other units 
such as pound-moles of product formed or gas pressure increase over a given 
amount of material for any volume and temperature. It is hoped that in this 
way the data obtained will be most useful to the largest number of people. 


Experimental Procedures 


Materials Handling 

Two simple gas tight propellant handling systems were built and 
tested for the purpose of filling containers for irradiation experiments. 

One system was used for fuels and the other for oxidizers. Figures 8 and 9 
are schematics of the propellant handling systems. Because of the corrosive 
nature of FLOX, OF 2 , the other oxidizers studied on this program, all the 
fittings, gages, tubing, etc., of the oxidizer handling system are constructed 
of stainless steel. Monel, or nickel. The fuel system is of brass and stain- 
less steel construction. Special care was taken in the cleaning and assembly 
of the components of the two systems to ensure removal of all contaminants, 
such as grease, moisture, oxide scale, etc., that might react violently on 
contact with oxidizers or fuels. The oxidizer system was passivated with 
fluorine gas which had been passed through an HF absorption tower before in- 
troducing FLOX or OF 2 . 

Each handling system contains a gas-measuring reservoir and pressure 
gage assembly of known volume so that a known amount of a propellant can be 
condensed into a precooled sample vessel for subsequent radiation effects ex- 
periments. The amount of propellant condensed may be controlled by succes- 
sively filling the reservoir and transferring its contents to the sample vessel. 
After a desired percent fill is obtained, the sample vessel valve is closed and 
the liquid sample and coolant container is loaded into the container in which 
the radiation experiments will be conducted. The handling systems also include 
facilities for transferring and loading liquid propellants, such as N 2 O 4 
studied in Task 1 and Task 3 of this program. 

Hydrazine (N 2 H 4 ) is the only propellant not handled in these systems. 
Anhydrous hydrazine is transferred directly into the sample vessels using a 
calibrated syringe to measure the volume transferred. The transfer is done in 
a simple glovebox arrangement under a blanket of dry argon in order to prevent 
the hydrazine from picking up moisture. 



FIGURES FUEL LOADING-SYSTEM 
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An additional gas handling system was built and leak-checked and then 
used to obtain PVT data on any gas buildup in the irradiated samples and to ob- 
tain samples for mass spectrometric analysis of radio lytic products. 

The fuel materials used in this study were analyzed by mass spec- 
trometry prior to their use in radiation experiments. 


FLOX . The fluorine-oxygen mixture was obtained from Allied Chemical 
and certified to be 70 ± 2 percent by weight fluorine. Our analysis by mass 
showed it to be quite pure. There was nitrogen present which could well have 
entered when we filled the sample vessel. A trace of hydrogen was observed. 
Some of the mass 28 peak could be due to CO since a peak was observed at 
mass 12 . 


Oxygen Difluoride . This material was obtained from Allied Chemical. 
Several analyses have been made and the purity is about 99 percent. There is 
a small amount of CO and HP present but their total quantity does not exceed 
1 percent. Nothing else was found in these samples. 


Nitrogen Tetroxide . The liquid N 2 O 4 was obtained from Air Products 
and Chemicals, Inc., and was found to be at least 99 percent pure. There was 
evidence of a small amount of NO and perhaps a trace of water. Nothing else 
was found and total impurity content was less than 1 percent. 


Diborane . The manufacturer, Callery Chemicals, claims the material's 
purity to be in excess of 98 percent and our analysis confirmed this. There is 
considerable nitrogen present which is to be expected, since nitrogen is gen- 
erally used to pressurize the tanks. A small amount of hydrogen was also 
found and evidence of trace quantities of higher molecular weight boranes. 


Hydrazine . Anhydrous hydrazine was obtained from Olin Matheson and 
found to haye purity >99 percent. Significant amounts of argon were found as 
a result of our procedure of filling the sample vessels under a dry argon 
blanket. Water was present as a trace constituent as was ammonia. There was 
a small amount of mass 28 which could be either N 2 or CO. 


Propane . This material was obtained from Air Products and found to 
be in excess of 99 percent pure. Only trace amounts of methane and ethane 
were found by mass analysis. 
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Sample Vessels 

The sample vessels for containing the condensed propellants for 
Task 1 and Task 3 radiation experiments are all of the same simple basic de- 
sign shown in Figure 10 but several materials were used. Some were con- 
structed using a Hoke 2HS10-305SS sample cylinder, a 7-inch-long 1/4-inch 
304SS high pressure tube extension, and a Hoke 1212G4Y valve for closure. 

Other vessels of the same type were made of Ti-6Al-4V alloy for use with 
hydrazine and N 2 O 4 in particular. Further, some experiments were conducted 
with hydrazine and propane in 1100 F aluminum vessels of the same design 
and dimensions as shoxm in the figure. 

Containers of 6063-T6 aluminum pipe construction in which the sample 
vessels were placed for carrying out radiation experiments at subambient and 
ambient temperatures in the Battelle gamma Co-60 radiation facility (swimming 
pool t 3 T)e) were built as shown schematically in Figure 11. A simple automatic 
liquid-nitrogen fill control system was constructed for Task 1 and Task 3 
radiation experiments conducted at liquid nitrogen temperature and extending 
over considerable periods of time. In these cases a Dewar was placed in the 
containers to hold the irradiation vessels in a liquid nitrogen bath. The 
containers were connected to the liquid nitrogen supply through rubber tubing 
insulated by two 1/2 -inch-thick lengths of Armflex insulating hose. This 
provided needed flexibility for inserting and removing the containers with 
sample vessels from the radiation source. Two thermistors were placed in 
the containers to sense liquid nitrogen levels. The lower thermistor would 
automatically start the fill cycle if the nitrogen fell below a predetermined 
level. The upper thermistor stopped nitrogen flow when the Dewar was filled. 
Flow was controlled by Asco solenoid valves. A timer was built into the cir- 
cuit to periodically add nitrogen to the Dewar even though the level had not 
fallen below that of the lower thermistor. 


Cobalt-60 Sources and Dosimetry 


Five Co-60 radiation sources were assembled and dosimetry measurements 
performed on each. The standard ferrous 'sulfate dosimeter was used. The dosim- 
etry measurements were carried out in a physical setup identical to that which 
was used in the Task 1 and Task 3 radiation experiments. The dosimeter was 
placed within a steel vessel of the same shape and wall thickness as the sample 
containers and this was placed in the aluminum container for irradiation. 

Thus, the dose rate seen by the dosimeter is the same as that which was seen 
by the sample. The dose rates (rads/hour) measured for the five Co-60 assem- 
blies were 6.0 x 10^, 8.0 x 10^, 5.8 x 10^, 5.0 x 10^, and 5.7 x 10^. The 
latter three sources were used for the long-term Task 3 experiments. Correc- 
tions were made for source decay as the program progressed. 


Cobalt-60 gamma radiation was used throughout the experiments since 
the radiation induced decomposition is, in general, independent of the nature 
of radiation and a function only of the amount of energy deposited in the fuel 
materials. Also, in an on-board situation, the major radiation to which a 
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fuel could be exposed may well be bremsstrahlung from the RTG. This is elec- 
tromagnetic radiation of the same nature as gamma rays. Thus, the results 
generated here should be generally applicable regardless of the nature of 
RTG in an actual flight situation. 


EXPERIMENTAL RESULTS 


Diborane (B 2 H 5 ) 

A series of experiments were performed with diborane in stainless 
steel vessels. All irradiations were at -78.5 C and a dose rate of 5.87 x 
10^ rad/hr. In most cases the sample contained 1.31 gm of B 2 Hg liquid, or 
about 3 ml, which is about a 20 percent fill of the vessel. The vessels 
were filled using the system previously described. 

The results of these experiments are given in Figure 12, where the 
moles of gas produced as a function of radiation dose is shown. The G-value 
for gas formation is G(H2) = 6.2 molecules/100 ev. This is equivalent to a 
gas buildup of about 0.14 ml/megarad. Assuming about one-half mole of gas 
is formed for each mole of diborane decomposed, the apparent rate of decom- 
position is about 0.04 percent per megarad. 

( 2 ) 

These results are consistent with those of Cornelius, et al., 
who found hydrogen to be the only gaseous product formed. The overall rate 
of decomposition may be somewhat higher in our experiments, but this is dif- 
ficult to ascertain accurately and probably not significant. The removal of 
a hydrogen atom from diborane would form an active species which would not 
be expected to remain long in the liquid without reaction. The products 
would be higher molecular weight borane compounds which will remain condensed 
under the conditions of these experiments. All such compounds might not be 
completely soluble but their rate of formation is probably not higher than 
that of hydrogen and hence their concentration is very low in these experiments. 

To further check this point, additional experiments with liquid di- 
borane were conducted in which mass spectrometric analysis was performed on 
both irradiated and unirradiated samples. As previously noted, the starting 
material is quite pure. In the vapor over the unirradiated liquid, BgHg 
and a very small amount of H 2 were the only species observed. There were 
traces of B and BH which probably arise from electron impact in the ionization 
chamber of the spectrometer. The sample was then completely vaporized and an 
aliquot taken. Again, no species other than B 2 H 5 and a small amount of H 2 
were observed. The same procedure was followed with an irradiated sample. 

The noncondensible gaseous product observed was H 2 . The irradiated sample 
was then completely vaporized and an aliquot analyzed. This showed com- 
pounds to be present which normally remain in the liquid when P-V-T measure- 
ments are made. These are the products expected from the dimerization of 
BgHg fragments resulting from irradiation if, as expected, one of the major 
initial effects is removal of a hydrogen atom from diborane. That is: 
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Other higher weight boron compounds would be expected to be formed in lesser 
amounts and this is also observed. The analysis showed a small amount of B 3 
compound and a trace of B 3 . 


To simulate storage under low dose rate conditions an experiment was 
carried out involving irradiation for a period of about 6 months (4376 hours) 
at the same temperature and in the same type of container as used in the higher 
dose rate experiments. In these long-term experiments the average dose rate 
was approximately 4.8 x 10^ rad/hr and the total exposure reached 22.2 megarad. 
The results of this experiment are also shown on Figure 12 where the result 
can be directly compared with those obtained at higher dose rates. It will be 
seen from the slopes of the lines that the apparent rate of radiation induced 
decomposition is not greatly different in the two cases and probably not sig- 
nificant, being 6.8 vs 6 . 2 . 


The curve for the low dose rate experiment does not go through the 
origin. The reason for this is not known. All sample preparation techniques 
were the same for the two sets of experiments. It is possible that the diborane 
warmed up during the 600 hours to reach the first data point. We do not know 
that this happened, but it is a possible explanation since thermal decomposi- 
tion increases rapidly as the temperature approaches room temperature. Control 
samples were prepared and stored under identical conditions, with the exception 
of radiation, and these showed a gas formation rate that was only about 10 per- 
cent that of the irradiated samples. However, the significant point is that 
the slopes of the two lines are essentially the same since this is a measure 
of the decomposition rate as a function of absorbed energy. Thus there does 
not appear to be any effect of dose rate over the range investigated in these 
experiments . 


An additional factor came to light too late in the program to be pur- 
sued in the data of Faust, Gould, and Kolman as reported in U.S. Patent 
3,001,920.(3) Although the technical data revealed is limited, it suggests 
that a higher conversion of diborane to products was observed by them than 
observed in our experiments over the same range of doses. However, apparently 
all of their experiments were done at temperatures from 0 to 40 C or possibly 
higher, whereas ours were at -78 C. This suggests there may be a significant 
temperature effect on the rate of conversion of diborane to products in a 
radiation field. This point should be investigated in any future work with 
this fuel. 
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Propane (C 3 H 3 ) 

Propane was chosen as the LPG fuel material of interest. This choice 
was made after discussions involving several NASA scientists as explained in a 
letter to the project technical manager. May 6 , 1969. Briefly, although there 
is considerable interest in methane as LPG fuel, there is also considerable de- 
tailed radiation effects data on methane. Propane, on the other hand, is also 
of interest as a fuel but there is, to our knowledge, no published radiation 
effects data on the liquid. Thus, it appears that, for this program, propane 
is the material which should be investigated. 

In summary, the most recent and complete work on liquid methane is 
that of H. A. Gillis.^*^^ He worked with the liquid at 112 K (-116 C) and used 
Co^*^ for the source of gamma radiation. The irradiations were carried to 
about 0.12 percent methane decomposition. The loss of methane per unit dose 
was given as G(-CH 4 ) = 6 molecules/100 ev. This is approximately 0.01 percent 
decomposition per megarad. 

The major products were hydrogen (H 2 ) and ethane (C 2 H 6 ) which to- 
gether account for about 95 percent of the observed products. Propane (C 3 H 8 ) 
was about 4 percent of the product and all other higher hydrocarbon products 
accounted for less than 1 percent of the total. The rate of their formation 
per unit of absorbed radiation is given by their G-values, listed below. 


G(H2) 

3.15 

G(C2Hg) 

2.21 

G(C3H8) 

0.23 

G (butanes) 

0.037 

G(all others) 

0.012 


Other radiation studies of methane have dealt primarily with the gas 
and have covered a range of temperatures in some detail up to 150 C, e.g.. Bone 
and Firestone , and Sieck and Johnson(6). Several other pertinent references 
are given in our report under NAS7-577, References 2-6. 


Experiments with liquid propane at high and low dose rates were 
initiated using the same general filling and handling procedures described 
previously. The irradiation vessels were filled in the usual manner using 
Ti-6A1-4V capsules of approximately 15 cc volume. All irradiations were at 
dry ice temperature. The short-term samples utilized a dose rate of about 
5.70 X 10^ rad/hr and the long-term experiments were at a dose rate of about 
4.80 X 10'^ rad/hr. At various exposures the buildup of gaseous products 
was determined by P-V-T relationships. Additional experiments with propane 
at the higher dose rate were done in 1100 F aluminum vessels to see if there 
were any obvious effects attributable to container material. 


The results of the experiments are shown in Figures 13 and 14. In 
Figure 13 the sample was irradiated at -78 C and measured at the same tempera- 
ture. In Figure 14, results shown are for samples irradiated as before but 
measured with the sample frozen at -196 C. At the lower temperature the most 
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likely decomposition products expected to remain volatile are H 2 and CH 4 . At 
the higher temperature, however, larger molecular species such as ethyl or 
butyl compounds would also be expected in the gas phase and would contribute 
to the total measured decomposition products. It is seen that apparent prod- 
uct yield as indicated by the G-value is 6.4 and 5.0 measured at -78 and -196, 
respectively. 

From the figures it will also be noted that the yield per unit 
weight of irradiated sample is somewhat smaller when the vessel is more nearly 
full of liquid. This is true at both measuring temperatures and the ratio of 
the 20 percent filled (1.76 gm) to 40 percent filled (3.96 gm) is about the 
same in both cases. The explanation of this result is not known. The results 
suggest that the propane vapor is more susceptible to radiation induced decom- 
position than the liquid. However, much more work will have to be done to 
determine if this is indeed the case. 

Samples of the gas above liquid propane at -78 C were analyzed by 
mass spectrometry. The results showed the major gaseous products to be un- 
saturated C 2 compounds, acetylene, and ethylene. Methane was also always 
present in significant quantities. The surprising result in these analyses 
was the general absence of hydrogen. From the amount of unsaturated compounds 
present there should have been easily discernible amounts of hydrogen. One 
might suspect that the hydrogen reacted with the titanium vessel in which the 
irradiation was performed. Such a reaction is less likely with an aluminum 
vessel, however, but there was no difference in the rate of product buildup 
when aluminum vessels were used. This can be seen in Figures 13 and 14 
where part of the data points are open circles and others are filled circles. 
The open circles are data obtained in titanium and the filled circles are 
data from aluminum vessels. There are no significant differences in the two 
cases. Also, an experiment was done wherein a liquid propane sample was 
irradiated to 13.5 Mrad in a titanium vessel and the pressure followed for 
98 hours after irradiation. No change was noted in the postirradiation period. 

Low dose rate experiments were also conducted with liquid propane 
in the titanium vessels for a period of almost 4 months at an average dose 
rate of 4.8q x 10^ rad/hr. These results are summarized in Figure 13 along 
with the data from higher dose rate experiments so that they can be compared 
directly. For the case of a 20 percent filled vessel the low dose rate data 
are indicated by crosses whereas the high dose rate data are shown by circles. 
The two sets of data fall on the same line indicating that the measured de- 
composition is function of the absorbed radiation dose but not the dose rate 
in this range. 


Hydrazine (N 2 H^) 

Preliminary radiation experiments with hydrazine in stainless steel 
vessels showed that considerable decomposition and gas buildup occurs. An- 
hydrous hydrazine ( 2.8 ml) was loaded under an argon atmosphere in a glove 
box into a 14.0-ml volume vessel. The sample was evacuated at -78.5 C through 
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several freeze -pump -thaw cycles to remove the argon and then irradiated at a 
dose rate of 6.0 x 10^ rads/hour at room temperature for 118 hours. As cal- 
culated from postirradiation P-V-T data, a gas pressure buildup of 44 psia and 
a 1.6 percent decomposition of the hydrazine occurred at the conditions of the 
experiments. Mass spectrometric analysis of the sample at room temperature 
showed considerable amounts of H 2 , N 2 , and NH 3 with an estimated molar ratio 
of 6 . 8:4. 3:1, respectively. The sample was then cooled to liquid nitrogen 
temperature, pumped on to remove the N 2 and H 2 , and allowed to warm to room 
temperature. Relatively large amounts of ammonia (NH 3 ), which had been con- 
densed while most hydrogen and nitrogen were removed, and trace amounts of H 2 , 
N 2 , and CO 2 were observed along with the hydrazine. No other radiolytic prod- 
ucts were observed. Preceding the mass spectrometric analysis, vapor pressure 
measurements of the sample obtained at liquid nitrogen. Dry Ice, and room 
temperature tend to support the mass spectrometric data. The starting hydrazine 
sample had a vapor pressure of about 3 psia (150 mm of Hg) at room temperature 
and only trace amounts of nitrogen and CO 2 were observed by mass spectrometric 
analysis. Further, preliminary experiments with hydrazine in stainless steel 
vessels indicated significant decomposition which seemed to be a function of 
the surface. Lucien and Pinns^^) had observed a surface effect when irradia- 
tions were performed in Pyrex vessels. Also, the apparent decomposition in 
the vapor was much greater than in the liquid at corresponding doses. 

A series of experiments was established for irradiating anhydrous 
hydrazine in Ti-6A1-4V vessels filled as described previously. All radiation 
exposures were at 22 C and a dose rate of 6 .I 5 x 10^ rad/hr. Several dif- 
ferent sample volumes were used so that the liquid-vapor ratio was varied. 
Hydrazine vapor was the only gas present over the liquid during irradiation. 

The results of these experiments are shown in Figure 15. The range 
of volumes of liquid samples was from 3 to 10 cc, which corresponds to 20 to 
66 percent of the total vessel volume. The figure shows that formation of 
gaseous products not condensible at -78.5 C, i.e., N 2 and H 2 , is unaffected 
by a change in ullage. The results are plotted as moles of gas formed per 
gram of sample as a function of dose to show that gas formation per unit of 
hydrazine is constant at the various gas-liquid ratios. The pressure build- 
up over the larger samples was of course greater. A larger mass of sample 
was being irradiated and the ullage was smaller leading to higher pressure. 

The results shown in Figure 15 are reasonably consistant with what 
other data are available. The rate of gas formation at STP is about 0.1 ml/ 

Mrad per gram of liquid hydrazine. In working with monomethyl hydrazine, 
which would be expected to give similar results, Shelberg(^) and Plank(5) 
both measured about 0.2 ml/Mrad. Lucien and Pinns give values between 0.1 
and 0.7 with most falling around 0.2 ml/Mirad. These values correspond to 
decomposition rates of the order of 0.01 percent per megarad. 

Based on volume of gas formed, our data give G = 5.8 (molecules of 
gas formed per 100 ev absorbed) . Based on the gas volume data of Lucien and 
Pinns we calculate G — 7 for their experiments with liquid hydrazine 
irradiated in glass ampoules by X-rays. The nitrogen-hydrogen ratios quoted 
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by them fall between 0.5 and 1.5 with most being 1.0 while our data give a 
value of about 1.5. 

Although this agreement is reasonably good it does raise some ques- 
tions. Lucien and Finns quote G-values for hydrazine decomposition of the 
order of 2 x 10^, which is difficult to reconcile with their gas release data. 
If hydrazine decomposes completely to nitrogen and hydrogen one would expect 
the gas ratio to be 0.5 and the G-value for decomposition to be one-third of 
that for gas formation, i.e., 


N2H4 “♦ N2 + 2H2 

However, if hydrazine predominately forms ammonia on decomposition the nitrogen- 
hydrogen ratio would be expected to be about 1 and the G-value for hydrazine 
decomposition would equal that for the formation of noncondensible gases, i.e., 

2N2H4 2NH3 + N2 + H2 

Our analyses have not shown large amounts of ammonia, nor apparently did those 
of Lucien and Finns. This matter needs to be investigated in further detail. 

An additional series of experiments was initiated with liquid hydra- 
zine in aluminum vessels for comparison with previous results obtained with 
liquid hydrazine in titanium alloy vessels. The initial results had indicated 
some surface-related decomposition in stainless steel vessels but not in the 
titanium ones. The later experiments were conducted using the same procedures 
throughout with only the vessel material being different. All experiments 
were at 5.53 ^ rad/hr and 22 C. 

The results of these experiments using aluminum vessels are shown in 
Figure 16. The volume of liquid samples used was either 3 or 7 cc, which is 
equivalent to about 20 or 45 percent of the vessel volume. The figure shows 
the moles of gaseous products not condensible at -78.5 C, i.e., N 2 and H 2 > as 
a function of radiation exposure. The products are given per unit weight of 
sample so that the two gas-liquid ratios for the different fillings can be 
directly compared. It is seen that decomposition, as determined by post- 
irradiation F-V-T measurements, is unaffected by a change in ullage. The 
G-value for gas formation is found to be 5.75 whereas the previously deter- 
mined value in titanium vessels was 5.8q. There is, therefore, no apparent 
difference between the two sets of data or the response to irradiation in the 
two types of vessels. 

Low dose rate experiments were conducted at an average of about 
5.6 X 10^ rad/hr for approximately 189 days (4544 hours) and a total exposure 
of 24.9 Mrad in the titanium alloy vessels. These data are summarized in 
Figure 17. The rate of product formation is equivalent to a value of G = 7.7 
which is significantly higher than in the high dose rate experiments in 
either type of vessel. Suitably stored control samples were maintained over 
the entire period of irradiation and the maximum increase in pressure in these 
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was less than 10 percent of the corresponding irradiated sample. The results 
were corrected for this nonirradiation pressure rise but it is apparent that 
a factor of this nature cannot account for the difference between the high and 
low dose rate experiments. The reason for the difference is not knovm and it 
does appear real. It is a point that should be considered in more detail in 
future work. 


FLOX (02 + F2) 

Experiments with liquid FLOX were conducted in stainless steel vessels 
at liquid nitrogen temperature at a dose rate of x 10^ rad/hr for a total 

exposure of 22.6 Mrad. The samples were prepared using the oxidizer handling 
facility and techniques previously described. The long-term tests with this 
oxidizer were allowed to progress for about 5000 hours. These samples received 
approximately 25 Mrad exposure. In none of these experiments have we seen any 
effect attributable to irradiation. This is not an unexpected result. The 
materials involved, liquid oxygen and fluorine, are stable compounds and any 
reaction initiated by the radiation would probably be small. The products of 
such reaction, perhaps 0 F 2 j would be present at low concentration, soluble in 
the liquid, and contribute very little to the vapor pressure of the sample. 

Due to the lack of observable effects, these experiments were carried no 
further. 


Oxygen Difluoride (OF 2 ) 


All short-term experiments with liquid oxygen difluoride were per- 
formed at a dose rate of 6.0g x 10^ rad/hr in stainless steel vessels and at a 
temperature of -196 C. Samples were loaded into the vessels using the pre- 
viously described handling system and pressure over the liquid was determined. 
Each vessel contained 5.0 gm of OF 2 liquid which filled about 20 percent of 
the available volume. After the appropriate irradiation, the pressure was 
again determined at -196 C and the moles of gaseous products formed were as- 
certained from P-V-T relations. 

The results of these experiments are shown in Figure 18 where the 
moles of gaseous decomposition products, O 2 and F 2 > are shoxm as a function of 
radiation dose. It will be noted that in this case the yield of decomposition 
products is not linear with dose as were the fuels. The fact that the yield 
is not constant but decreasing probably indicates a secondary reaction which 
is removing the gaseous products at a slow but definite rate. One such re- 
action which immediately suggests itself is the reaction of O 2 and F 2 gas to 
reform OF 2 liquid, i.e.. 


1/2 02(g) + F2(g) 0F2(1) 

It is also possible that some O 2 F 2 could be formed which would also condense 
at the temperature of these experiments and lead to an apparent decrease in 
the rate of decomposition of OF 2 • 
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Since the rate of formation of decomposition products, and hence 
the G-value, decreases with increasing dose, the value at low doses is most 
important since that gives the maximum value and indicates the maximum rate 
of radiation induced decomposition. The initial value is G(gas) = 0.30 
molecules/100 ev. A G-value this low indicates liquid OF 2 is rather radia- 
tion resistant. The data further indicate a maximum gas buildup rate of 
about 6.9 X 10“^ ml/megrad. Assuming one and a half molecules of gaseous 
product for each molecule of liquid OF 2 decomposed, G-( 0 F 2 ) = 0.20 molecules/ 
100 ev. This is equivalent to about 1 x 10"^ percent decomposition per 
megarad. 


Long-term experiments were extended for a period of almost nine 
months. However, after about five months we began to experience difficulty 
with leaking valves on the sample vessels so that the data after that were 
not reliable. As can be seen from Figure 18, the rate of buildup appears 
more linear at 5.6 x 10^ rad/hr than at the higher dose rate. However, the 
amount of decomposition is always very small in all experiments and the dif- 
ferences in the two cases are probably not significant and using an overall 
value of G(gas) = 0.17 for the whole range of dose rates is justified. 


Nitrogen Tetroxide (N 2 O 4 ) 


Experiments with liquid nitrogen tetroxide were conducted by the same 
general procedures outlined previously. Stainless steel vessels of approxi- 
mately 15 cc were used. Irradiations by Co-60 were performed at a dose rate 
of 5.76 X 10^ rad/hr and a temperature of 21 C. Sample vessels were evacuated 
prior to filling with N 2 O 4 . After filling, the liquid N 2 O 4 was frozen and any 
noncondensible gases removed by pumping. Noncondensible gases at liquid N 2 
temperature formed by decomposition of N 2 O 4 were determined as a function of 
time and radiation dose by P-V-T relationships. 

The results of experiments at the dose rate of 5.7g x 10^ rad/hr for 
both the uninhibited and inhibited forms of N 2 O 4 are shown in Figures 19 
through 22 as moles of noncondensible gas formed per gram of liquid N 2 O 4 as a 
function of radiation dose. Several points are immediately obvious from these 
figures; the inhibited form shows more apparent radiation decomposition at a 
given radiation exposure; the scatter in the data, particularly at the lower 
filling (3-cc sample), is greater than in many of our other experiments; and 
the decomposition rate at a fixed dose rate appears to be a function of the 
extent to which the sample vessel was filled. A 3-cc sample fills the sample 
volume approximately 20 percent while a 10-cc sample fills it about 67 percent. 

These results suggest that the vapor is more susceptible than the 
liquid to decomposition and that there is some interaction with the walls of 
the vessel. Several other observations are consistent with this. The amount 
of noncondensible gas formed, largely N 2 > in unirradiated control samples was 
also a function of the filling, being between about 1 and 20 percent of the 
radiation induced decomposition product for the 3 -cc samples, and 2 percent 
or less in those samples with 7 o r 10 cc of liquid N 2 O 4 . The magnitude of 
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these numbers are only significant for the particular dose rate used, of 
course, but the observed differences between samples again suggests an inter- 
action between the walls and nitrogen tetroxide vapor. It was also noted 
that the rate of the nonradiation decomposition, as indicated by buildup of 
gaseous products, was greater immediately after filling than at later times. 

In several experiments all of the gases were examined rather than 
just the noncondensible species. The major component, excepting the N 2 O 4 , was 
NO in both the irradiated and nonirradiated cases and its rate of formation 
appeared about the same in both cases at the dose rate used, 5.?5 x 10 ^ rad/hr. 
Nitric oxide is an expected radiation decomposition product and a reasonable 
product of an oxidation reaction of NO 2 with the walls of the vessel. 

It will be noted that the relative rates of radiation and nonradia- 
tion reactions is a direct function of the dose rate. At dose rates lower 
than those used in these experiments the effect of the radiation could be a 
small fraction of the overall reaction at least when storage vessels were not 
full of liquid. 


The only significant noncondensible gas found on mass spectrometric 
analysis was nitrogen. The G-value for noncondensible gas formation then is 
the G-value for N 2 formation. For the uninhibited form G(N 2 ) = 6.1 x 10“^ 
molecules /100 ev which can be compared to the results reported by Castorina^^®) 
of 5.2 X 10“^. As seen in Figure 19, the G-value for the inhibited form is 
higher, being G(N 2 ) = 18.1 x 10"^ for the 3-cc samples. A summary of the 
G-values for noncondensible products at the various loadings is given below. 


3 cc 
7 cc 
10 cc 


N 2 O 4 Inhibited 

18.1 X 10"2 
5.4 X 10"2 
3 . 6 X 10”2 


N 2 O 4 Uninhibited 

6.1 X 10"2 
2.8 X 10-2 

2.2 X 10'2 


The only other product found in addition to NO and N 2 was nitrous 
oxide (N 2 O) which is formed at a rate slightly greater than that of nitrogen 
in the case of the uninhibited nitrogen tetroxide and slightly less in the 
case of the inhibited form. This again is consistent with the results of 
Castorina, who reported the formation of N 2 O to be about half of that of N 2 • 


The lack of oxygen among the products is probably due to the NO 
formed, even in the absence of radiation. Any oxygen formed by decomposition 
of NO 2 by irradiation would react rapidly with NO to reform NO 2 or react with 
the wall which could act as a sink for the oxygen. 


In general, the differences between the short-term and long-term 
experiments do not seem to be very large as can be seen from inspection of 
Figure 22. This indicates that the radiation decomposition in these cases 
is probably not strongly dependent on dose rate and the extent of decomposi- 
tion will always be directly proportional to the radiation energy absorbed by 
the propellant material. The G-value corresponding to the line in the figure 
is 21 X 10"2 for a 3-cc sample compared to 18 x 10“2 for a 3-cc sample of the 
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inhibited N 2 O 4 at high dose rate. The G"value for the uninhibited material at 
the low dose rate is somewhat higher, however, than the corresponding high 
dose rate case. Such differences as are seen may occur because thermal decom- 
position becomes a much larger factor at the lower dose rates. Although un- 
irradiated "control" experiments have been carried out for each material, 
uncertainty from this source becomes more significant in the low dose rate 
cases . 


The apparent rate of the nonradiation induced gas product formation 
decreases with increasing exposure time. On the other hand, the radiation 
induced decomposition occurs at an essentially constant rate proportional to 
the dose rate. At the dose rate of these experiments, approximately 5.6 x 10^ 
rad/hr, the thermal reaction predominated during the first 6 weeks of exposure. 
After that the rate of the thermal decomposition had fallen to the point where 
the radiation induced decomposition becomes the dominant factor in the buildup 
of noncondensible gases. With continued irradiation, the buildup of gaseous 
products per unit absorbed dose is about the same as observed in the high dose 
rate case, the inhibited N 2 O 4 shows a somewhat higher rate of radiation decom- 
position than the uninhibited N 2 O 4 . Because decomposition is slow, the pressure 
buildup was only a few torr per month and the thermal contribution during the 
first couple of months was comparatively large. A small error in this correc- 
tion could possibly account for the curve, as shown in Figure 22, not going 
through the origin as expected. 


Summary of Radiation Effects Data 

The most striking result is the obviously much larger sensitivity to 
radiation induced decomposition of the fuels, as a group, compared to the 
oxidizers. Although molecular oxidizers are more susceptible than elemental 
ones, e.g., FLOX or liquid oxygen, they are still much more stable in a radia- 
tion field than the fuels examined in this program. This has obvious implica- 
tions in regard to the positioning of fuel and oxidizer storage containers with 
respect to nuclear power sources and any shielding incorporated in the system. 
Although there are some differences among the materials studied with regard to 
container material and dose rate, as discussed in the previous sections, a 
rough comparison among these materials is shown in Figure 23. Although such a 
comparison cannot be exact, it is nonetheless instructive. 


Using data on isotope sources previously reported (Contract NAS7-577), 
off-gas volumes formed in various propellants by plutonium, curium, and stron- 
tium heat sources are estimated in the following tabulation. For purposes of 
the estimates: the sources are estimated to be 200,000 watt (th) at 1 foot for 
°Pu 02 and ^^Cm 203 , which would be equivalent to 10 kw(e) at 5 percent con- 
version. The corresponding estimated dose is 1.3 x 10^ rad, ^^^Pu 02 and 
1.1 X 10^ rad, Cm 203 . A 2000-watt (th) (100 w electrical) ^®SrTi 03 source 
at 1 meter is assumed to give a dose of about 3 x 10^ rad. 



RELATIVE GAS VOLUME FORMED PER UNIT PROPELLANT PER UNIT DOSE 
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Estimated 

Volume, ml. 

of Off-Gas 


Per Gram 

of Propellant 

Per Year 

Propellant 

^^^Pu02 

2‘^'^Cm203 

^°SrTi 03 

B 2 H 6 

0.03 

0.2 

6 

C3% 

0.02 

0.1 

5 

N2H4 

0.02 

0.1 

5 

FLOX 

0 

0 

0 

OF 2 

0.0003 

0.0002 

0.06 

N 2 O 4 (uninhibited) 

0.0003 

0.0002 

0.06 

N 20 ^ (inhibited) 

0.0006 

0.0004 

0.1 


The radiation doses used in the above estimates assume no self- 
shielding or external shielding. Since self -shielding will be significant in 
all of these sources at the power levels assumed, the estimates can be treated 
as maxima. It appears that the a-sources would present no large problems with 
these propellant materials over periods of several years. However, the ^®Sr 
p-sources might present some pressure buildup problems, at least in the case of 
fuels. From the experiments completed to date, it appears that the pressure 
would continue to build up over the propellant linearly with time over mission 
times of 1 to 10 years. The rate per unit time will depend on the size, con- 
figuration, and shielding of the radioisotope source. If the source decays 
significantly during the mission, a corresponding decrease in the rate of 
pressure buildup would be expected. That is, total noncondensable gas re- 
mains proportional to total dose. The accumulation of gas can be a signifi- 
cant engineering problem even though the amount of liquid decomposed by the 
radiation is relatiyely small and, of itself, probably not significant. 
Determining how best to accommodate those gaseous products which are formed 
could require a major technological effort. 

Because there are indications of some interactions resulting from 
surface effects with particular propellant materials and specific container 
materials as well as some unexplained differences at low dose rates, addi- 
tional research probably will be needed on specific propellant systems as 
these become better defined for given missions. This appears to be particu- 
larly true for the hydrazine and diborane fuels and the nitrogen tetroxide 
oxidizer. Interactions with container materials may well lead to pressure 
buildup of a similar nature to that from radiation. For long missions where 
the radiation dose rate might be expected to be minimum, surface interactions 
with the container might well be the major source of gaseous products. Shield- 
ing of the radiation source would then not accomplish the desired reduction in 
gaseous product accumulation. Further, the rate of such nonradiation decom- 
position would be expected to be a much stronger function of temperature than 
is the radiation induced reaction. The influence of this factor might only 
become apparent at low dose rates. Thus, further research should be directed 
at the elucidation of these factors. Only a more thorough understanding of 
the mechanisms, both radiation and nonradiation, and their interaction will 
permit more reliable predictions of gas buildup to long periods of time. 
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APPENDIX A 


DEFINITION OF COMMONLY USED TERMS (NOMENCLATURE) 


Split Shield . A shield configuration consisting of two separate 
shields, one between the radiation source and the propellant tanks and one 
between the propellant tanks and the instrument package. 


Laminated Shield . A shield configuration consisting of multi- 
layered shield component (s) . 


Shield Segment . One of the shields in the split-shield 
configuration . 


Primary Direct Shield . The shield segment between the radiation 
source and the propellant tanks. 


Secondary Direct Shield . The shield segment between the propellant 
tanks and the instrument package. 


Scattering Structure . Any structure which may reflect a portion of 
radiation incident upon it to the propellant tanks or instrument package. 


Scatter Shield . A shield used for the purpose of reducing the 
radiation reaching a scattering structure. 
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APPENDIX B 


TRANSMISSION MATRIX DERIVATION AND USE 


A Laplace transformed transmitted flux function may be written; 
F(S) = H(S) • G(S) 


where 


F(S) is the Laplace transform of the transmitted flux function, 

G(S) is the Laplace transform of the incident flux function, and 

H(S) is the Laplace transform of the transmission function. 

Since it is not intended to determine the Laplace transformation of all source 
distribution which may occur, the simplest and most useful source function, a 
delta function in energy denoted 6(E), is chosen. The Laplace transform of 
the delta function is unity and, therefore, simplifies the analysis consider- 
ably since any source spectrum may be described as a continuum of discrete 
energy sources. Thus: 


F(S) = H(S) • 1 
or F(S) = H(S) . 

This implies that the inverse Laplace transforms of F(S) and H(S) are identi- 
cal. Therefore, if discrete excitation energy-group transmission data are 
available, it is possible to construct the transmission function for any 
arbitrary input spectrum by summing up the discrete excitation transmission 
functions . 

The advantage of using the transmission functions is that the group- 
to-group transfer matrix which is made up of point values from the transmission 
function is totally described for a given slab thickness. Thus, iterative sol- 
utions are not needed and significantly less computation time is anticipated. 
For laminated or multilayered shield configurations, this technique eliminates 
the uncertainty in buildup factor determination in the case of gamma radiation 
and also eliminates the necessity for use of fictitious removal cross-section 
information in the case of neutron radiation. 

The most important assumption inherent in the technique is that 
double reflections across material boundaries are insignificant* contributors 
to the final solution. 


* Ph.D. thesis, J. L. Ridihalgh, Iowa State University, 1968. 
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APPENDIX C 


LISTING OF THE OPTIMIZATION COMPUTER PROGRAM 
FOR THE PRIMARY DIRECT LAMINATED SHIELD 


A listing of the optimization computer program for the primary 
direct laminated shield appears in this appendix along with sample-problem 
input data and output. The computer program consists of seven routines. 
They are: 

(1) The main program which performs the iteration- 
control functions necessary to properly opti- 
mize the shield system 

(2) Subroutine Geom which sets up the geometry in- 
formation necessary for the program operations 

(3) Subroutine WATE which determines the weight of 
any shield system designed by the program 

(4) Subroutine Transf which reads radiation trans- 
mission factor information for the source and 
shield components 

(5) Subroutine Flux which computes the radiation 
transmitted through the shield to the pro- 
pellant tank 

(6) Subroutine Source which computes the source 
radiation arriving at the front face of the 
shield 

Subroutine Intrp which is used by Subroutines 
Source and Flux to perform the interpolations 
and/or extrapolations on the transmission 
factor data. 


( 7 ) 
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Sample Problem 

The sample problem is a two-group two-laminate shield optimiza- 
tion problem. The transmission factors used are purely fictitious since 
the problem was set up only to demonstrate the program operations. The 
order of input data for this problem is: 


(1) Diameter of source - 5 cm - DSRC 

(2) Length of source - 10 cm - XLSRC 

(3) Diameter of propellant tank - 30 cm - DTANK 

(4) Separation distance between propellant tank and the 
radiation source - 100 cm - SEPAR 

(5) Number of shield laminates - 2 - LAM 

(6) First guess of laminate thicknesses - 5.^ 5. - X(I) 

(7) Density of shield laminates - ,2., 1. - RH0(I) 

(8) Number of energy groups - 2 - N0G 

(9) Dose conversion factors - l.E-4, l.E-5 - FD(I) 

(10) Limiting dose rate - 1. - DRLIM 

(11) Source strengths/vol - l.E+10^ 0. - SRCF(I) 

(12) Number of thicknesses for which the transmission 
factors are inputed - 4 - NTHK 

(13) . Thicknesses for the shield laminates for which the 

transmission factors apply - 

First Laminate Second Laminate 


1 cm 

2 
8 

16 


1 cm 

2 
8 

16 


1, 1, 2, 2, 8, 8, 16^ 16 - TKTLM(I,J) 
where I = Laminate Numbers 1 and 2 
J - Thickness Numbers 1 to 4 
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(14) Transmission factors for the source material = 

First Thickness 


0.9 

0. 

0.1 

0.95 

Second 

Thickness 

0.8 

0. 

0.2 

0.89 

Third ' 

Thickness 

0.6 

0. 

0.4 

0.8 

Fourth 

Thickness 


0.3 
0.7 

where the matrix is 

Position 11 - Group 1 

Position 12 - Group 1 

Position 21 - Group 2 

Position 22 - Group 2 

- TSRC (I,J,K) 

where I = receiver group number 
J = incidence group number 
K = thickness number 

(15) Thicknesses for the source for which the transmission 
factors apply - 1, 2, 4, 8 - TKSR(I) 


0 . 

0.6 

11 12 

21 22 

transmission for Group 1 incidence 
transmission for Group 2 incidence 
transmission for Group 1 incidence 
transmission for Group 2 incidence 


where 1 = thickness number 
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(16) Transmission factors for the shield laminate materials 


Material for First Laminate Material for Second Laminate 

First Thickness 


0.5 

0 

0.9 

0 

0.49 

0.99 

0.1 

0.5 


Second 

Thickness 


0.25 

0 

0.8 

0 

0.72 

0.98' 

0.19 

0.25 


Third Thickness 


0.02 

0 

0.5 • 

0 

0.9 

0.9 

0.4 

0.02 


Fourth 

Thickness 


0.0005 

0 

0.25 

0 

0.8 

0.8 

0.2 

1 O.OOl 


The actual input and resultant program output are listed 
below. 


Input 


SDATA 

10.^ 30* j 1 00*^ 5* j 5. ^ 2« ^ 1 • 2^ • 000 W J • E- 5^ 1 . ^ 1 • E+ 1 0/ 

0* 1 1 2. # 2«^ ^ 1 6* J 1 6» • 9^ 0» ^ • B# 0. j . 2^ *89# 

. 6# 0. > • • 8j - 3^ 0»i . 7# ♦ 6^ 1 * j 2* J 4. ^ 8* J » 5^ 0* / • '^9^ • 99# » 9# 0« # 

. 1# • 5# • 25# 0» # • 72# •98# .8^0»#*19#»25# •02* 0»#«9#»9#«5#0»# 

• A.* » 02# .0005# 0.# «8# *8# ->25# 0.# .2# • 0005 


Output 


ITERATION 
1 5.1 661 
1 4. 1006 
15. 18 08 
■ 16.5302 
20. 5709 

laminate 

i 6. 4003 


MONITOR 

1 5. 1 661 
1 7.2341 
1 5. 1808 
1 3. 2241 


12. 7539 
THICKNEGSE5 
13. 1202 


1 064 . 5226 
1 1 45. 7562 
■1066. 146 
1025. 9523 
1233.671 

WEIGHT 

1025.9523 




10 COMriON RH0C2) 

20 COi-.KON XC2)^ Fr C2)> SRCF < 2 ) ^ TL AMC 2> 2# 2> 4)# TSKC< 2^ 2:» 4)> T!<TLMC2^ / 
30 -rlKSFCC RATIO (2)^ SLOP^ PuMIN# DSRC> DTANK# LAWj XLSRC^ SEPAR^ WF.I GH'f^ 
40 -J-NOG ■ • ■ . 

50 1 

60 CALL GEOM 

65 LEI GHT= 1 00000. ' • 

80. REAI^NOG^ CFKD^ I = 1,N0G)> DRLIM 
9 0 READ# CSRCFU)^ I-1,N0G) 
no CALL TRAN SFCNTHX) 
in PRINT#" ITERATION MONITOR" 

130 12 CALL FLUX(NTHK# DOSE) 

MO IF(DRLIM-D0SE)9 
1 50 IFCDRLIM-. 1 '^KL IM- DO SE) 10 
151 IFC D0SE-DflLIM)20 
155 YY»0* 

160 DO n I=1#LA.M 
165 YY=YY+X(I) 

170 11 xa )=XC I )=!=. 9 

175 IFCYY-.D18 

190 60 TO 12 

191 20 DO 21 I=1#LAM 

192 21 X(I)=X<I)=!=.9S 

193 GO TO 12 

200 9 DO 13 I = 1#LAM 

210 13 Xn)=X(I)-M. 12 

220 GO TO 12 
230 10 XX=0* 

240 to 14 1=1# LAM 
250 14 XX=KX+XCI) 

260 DO 15 1 = 1# LAM 

270 15 RATIOC I )=X< I )/XX 

280 CALL WATECWT) 

290 IF<2“IJ) 18 

295 PRINT# CXCI )# 1= 1#LAM)# WT 

300 I F C 1 . 0 1 ■!•• EE 1 GH T - WT ) 1 9 

310 WEIGHT=WT 

320 IFC 1-IJ) 1 7 

330 IFC.99-RATI0C 1)) IS 

340 IJ= 1 

350 RATIOC 1 )=.9ARATI0< 1 ) 

360 RATI0(2)=1. -RATIOC 1 ) 

370 XC 1 ) = RATI0C 1 )=i‘XX 
380 XC2)=XX“XC1) 

390 IFCXC 1 )-. 01*XX) 13 

400 GO TO 12 

^10 19 IFC IJ-2) 1 7 

420 GO TO 16 

430 17 IJ=2 

435 WEIGHT=WT 

4 40 XC 1 ) = RATIOC 1 ):^<XX/.9 
450 XC2)=XX-XC 1 ) 

^60 IFCXC2)- . 01-XX) 18 
470 GO TO 12 

480 1 6 XC 1 ) = RA7 IOC 1 )>i^XX=:^. 9 

490 XC2>=XX~XC i ) 

500 IJ=3 
510 60 TO 12 

520 13 PRINT#" LAMINATE THICKNESSES WEIGHT" 

5 3j0 PR I N T^^X C I ) # I =j_# L A M ) # J;^E I GH T 
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55 0^ S UFiPOUTINfe C-F.OM ■ 

S55 C6i'-n'.0i\' RHOCR) 

5 60 COMMON XC2)> F DC 2 ) ^ SRCF < 2 ) :» TLAM ( 2:- 2:. 4)7SKC(2^ 2^ /i)j TK7LM(2^ 4)^ 

5 65 +'1KSRCC 4) j> F^AT 1 0 C 2 ) ^ SLO F:M I N> DSRC> DTANK^ LAM^ XLSHC* SEPARj WEI GHT^ 
570 +NOG 

5 75 READ;. DSRC^ XLSRC^ DTANK;. SEPAR^LAM> <Xd )> I = 1>.LAM) 

6 1 5 READ^ ( RHO C I) d = 1 ^ LAM ) 

630 ' IFC DSRC-DTAN10 2 

635 SL0P=. 5=dDTANK-DSRC)/<SEPAR+. 5=^=DTANK) 

640 F.MIN=. 5RDSRC+SL0p::dEPAR 
6.fi5 GO TO 3 

650 2 SLOP-. 5=FCDTANK-DSRC)/(SEPAR+. 5H;DTANK + XLSRC) 

655 RMIN=. 5=!^DSRC+SL0P-i:XLSRC 
660 3 CONTINUE 

670 CALL WATECWT) ' 

67 6 WEIGHT=^WT 

680 RETURN 

700 SIJBR!^' I N E JvAT^C_d ) , 

7f0 MMErLSIoT ■A(2)dOL(2) 

715 COMMON RHO (2) 

720 COMMON X ( 2 > ;. FD< 2 ) ;. SRCF ( 2 ) > TL AM C 2^ 2^ 2^ TSRCC2> 2d)> TKTLMC2^ 

725 +TKSRCC 4)dATI0(2)dL0P^ RMIN> DSRC/ DTANK^LAM^ XLSRC^ SEPAR^ WEI GHT^ 
730 +N0G 

735 WT=0. , 

740 RMINl^RMIN 

7 45 DO 4 K=;1^LAM - . 

750 IFCK-2)5 

755 A<K) = ACK- 1 ) + £LOP>i:X(K) 

7 60 GO TO ■ 

765 5 IF(DSRC-DTANK)6 

7 70 DO 7 1=1 dAM 

775 7 KMIN 1 = RMIN-X( I )-tSLOP 

780 6 AC1)=RMIN1 

785 4 CONTINUE 

790 DO 8 K=1>LAM 

79 5 VOL(K) = C ACK)N.:’,<2)=!=X(K) + SL0P*C ACK)4:X(K)=:-=H:2+SL0P/3.*X(K)+A-3) 

800 8 WT=VOL<K):FRHOCK)+ WT 

80 5 RETU RN 

8 5~0~ S y B.ROLI T I Nj; TRA NS F(N THK) ^ 

8 55 C0MM0N“rH0(2) 

8 60 COMMON X(2)j FD C2) ^ SRCF ( 2 d TL AM ( 2^ 2> 2 j /id TSRCC 2d^ 4 ) dK TLM C 2» / ) ^ 
8 65 +TKSKC( s RAT 1 0 C 2 ) ^ SLO P^ RMIN, DSRCj DTANK;.LAM;. XL £RC>' SEPAR^ WEIGHT^ 
8 70 +N0G . 

875 READ^NTHK^ ( < TK TLM C I J ) , I = d L AM ) d= I > NTHK ) 

880 READ^ CCCTSRCCI^ I = d NO G) > J= d NO G) ^ K= d N THK ) 

885 READdTKSRC(I)d=dNTHl<) 

890 READ^ <C( (TLAMCd Jd<d)^ I = d NO G) ^ J= d NO G ) d= d L AM) d= d N THK > 

895 RETURN 
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9 10 SUBROUT K'j'E FLUXCi'nilK^ DOSF.) 

9 T 5 COl'WON Rj^b'Cg) 

920 COMMON XC 2 ) > F D< 2) , SRCFC 2 ) > TL AM < 2> 2^ 2> TSRC(2^ 2^ /i)^ TKTLMC2/ AO j 
9 25 + 7i( SF?C ( / ) ^ RAT 10 ( 2 ) ^ SLO PM IN^ DSRC^ D'l ANK> LAM, XL SRC, SFPAR, IvEI GHT, 
930 +NOG 

935 DIMENSION FLUC 2) , FLUX < 2 ) 

9 AO DO 3 1=1, NOG 

945 FLUX<I) = 0. ■ 

950 3 FLU(I) = 0. 

9 55 ' GO TO 8 
9 60 5 DO 2 K=1,LAM 

965 DO 1 1=1, NOG 
9 70 DO 1 J=I,NOG 

975 CALL IN7RP( I, J,K,XCK), TL,NTH:<) 

980 FLUXC J) = FLU( I)--:=TL+FLUXCJ) 

985 1 CONTINUE 

990 DO 2 1=1, NOG 
995 FLUC I ) = FLUXC I ) 

1000 FLUX(T)=:0. 

1005 2 CONTINUE 

1010 DOSE=0. 

1015 DO 6 1 = 1, NOG 

1020 6 DOSE=DOSE+FDC I )n--FLUC I ) 

102 5 rOSE=DOSE/C 12. 566*SEPAR-F*2) 

1030 RETURN 

1035 8 CALL SOURCEC FLU, NTHK ) 

1040 GO TO 5 • 

J 0_£5 P^.TUt<il ^ 

foso SUD~ROUTINE S O UR C EC FL U, NTHK ) 

1 055 COMMON RH0C2) ~ - - - . 

1060 COMMON XC2), FDC 2 ) , SRCF C 2 ) , TLAMC 2, 2, 2, 4), TSRCC2, 2, 4), 

1065 +TKTLMC2, 4), TKSRCC 4),RATI0(2), SLOP, RMIN, DSRC, DTANK, LAM, XLSRC, 
1070 +SEPAR, WEI GHT,NOG 

1075 DIMENSION FLUC 2) ^ FLUXC 2 ) , FRUC 2) 

1076 DO 1 1=1, NOG 

1077 1 FKU(I)=0. 

1030 DO 4 M=l, 10 

1085 XSRC=XLSRC/ iO.H=CM- 1 ) + XLSRC/20. 

1090 DO 3 K=1,N0G 

1095 3 FLUCK) = (SRCFCK)i!'. 785:i'DSRC’i-*2)=i'XLSRC/ 10. 

1100 DO 2 1=1, NOG 
1 105 DO 2 J=I,NOG 

1110 CALL INTRPCI, J,0,XSRC, T£,NTHK) 

1115 FLUXC J) = FLUC I)--KTS+FLUX(J) 

1 120 2 CONTINUE 

1125 DO 4 1=1, NOG 

1130 FRUC U = FLUXC I ) + FRU( I ) 

1135 4 FLUX(I)=0. 

1140 DO 7 1 = 1, NOG 
1142 7 FLUCI )=FRUC I ) 

T145 RETURN 
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1 1 50~ SUB ROUT I'M i MTRF ( I > K> X 1 > T> N ) 

U 5 5 CO FiMO N RHO (‘g ) 

1 160 COMMON X(2)^ FD(f>)^ 5RCFC2)^1LAM(2>{?^2^ /i) TSFCC 2^ 2> 4)>TKTLM(2# 4)^ 
1 1 65 4-7KSRCC 4)^ KATI0(2);» SLOPj RMIN> DSRC^ DTANK> L AM^ XLSRC^ SF.PAR> 

1170 +WEIGH7>N0G 

1175 IF(K”I)1 

1 IhFO DO 2 II = l:»N 

1135 2 IFCX1“TKTLMCK> II ) )3 

1190 Z = X1"TKTLMCKW'J) 

1195 H=TKTLM<K,N)-TK7LM<K^N- 1 ) 

2000 R^TLAMCI, J>K>N)-TLAiM<I> J^K>N- 1 ) 

2005 7-tLAM(I> JjK^N) + R/H^=Z 

2006 IFCDll 

2007 IFC1.-T)12 

2010 RETURN 

2011 11 T=0*i RETURN 

2012 12 T=l.i RETURN , 

2015 3 IF(II-2)4 

2020 Z=CXl-']KTLMCKi I I- 1 ) )/CTKTLM(K> 1 1 )-TKTLM(K> II-l ) ) 

2025 H=TLAMCIj J^K> II)-TLAM<I> J>K^ II- 1) 

2030 T=TLAM( I, 1 1- 1 ) + Ht-Z 

2035 RETURN 

2040 4 IF<J-I-1)9 

2050 T = X1/TK1LM(KM )^TLAMC I;. J>K> 1) 

2055 RETURN 

2056 9 T= 1 . -XI/TKTLMCKj 1 1 • -TLAMC I^ K> 1 > ) 

2057 RETURN 

2060 1 DO 5 II = 1^N7HI< 

2065 5 IFCXl-TKSRCCII ) 1 6 

2070 Z=X1-TKSRCCN) 

2075 H=TKSRCCN)-TK5RCCN- 1 ) 

2080 R=TSRC(I^ J>N)-TSRC<I^ J>N-1) 

2085 T=TSRC( I^ JjN)+R/H*Z 

2036 IFCTlll 

2087 IFC i.-T) 12 _ 

2090 RETURN 
2095 6 IF(II-2)7 

2100 Z=CX1-7KSRCCII-1)1/CTKSRC(II)-TKSRCCII-1>) 

2105 H=TSRCCn II)-TSRCU>J^ II-1> 

2110 T=TSRCCIjJ^ II-U+H'-:-Z 

2115 RETURN 

2120 7 IF(J-I-1)13 

2130 T = X1/7KSRCC D^i^TSRCCI^ J> 1) 

2135 RETURN 

2136 13 T=lo-Xl/TKSRCCl)>;:( l.-TSRC(I>JM) ) 

2137 RETURN 
2140 END 
2200 SUSE DATl 
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